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ZYAM and V2: Underestimating jet yields from dihadron azimuth correlations 
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Dihadron azimuth correlations can provide combinatoric access to jet structure in nuclear colli- 
sions. To isolate true jet yields a background must be subtracted, including a constant offset and 
a contribution from "elliptic flow" (azimuth quadrupole measured by i;|). The principle of "zero 
yield at minimum" (ZYAM) has been introduced to determine the constant offset. Independent 
measurements determine V2- This analysis demonstrates that the ZYAM concept is invalid (offset 
typically overestimated) and V2 is also overestimated by conventional measurements. Jet yields are 
thus substantially underestimated in more-central A-A collisions, and the "away-side" azimuth peak 
(back-to-back jet correlations) is strongly distorted, leading to incorrect inference of "Mach shocks." 

PACS numbers: 12.38.Qk, 13.87.Fh, 25.75.Ag, 25.75.Bh, 25.75.Ld, 25.75.Nq 



I. INTRODUCTION 

The nominal goal of the Relativistic Heavy Ion Collider 
(RHIC) is production of quark-gluon plasma (QGP), a 
color-deconfined, thermalized state of quarks and glu- 
ons Several characteristic indicators for QGP for- 
mation have been proposed, one being modification of 
parton scattering and fragmentation to hadron jets in 
more-central A-A collisions [2j]. Colored partons should 
be uniquely sensitive to the presence and properties of a 
color-deconfined medium. In A-A collisions jets are stud- 
ied combinatorially through dihadron correlations ana- 
lyzed and interpreted in_the context of parton energy 
loss in a dense medium 3]. Two diametrically-opposed 
pictures have emerged from jet correlation analysis. 

In the conventional picture the scattered-parton spec- 
trum extends down as far as 1 GeV, resulting in copious 
parton production near that energy (estimating ~1000 
gluons per unit rapidity) [1, Q . Most partons deposit a 
large fraction of their energy in a "medium" and are ther- 
malized Q , producing a large energy density with ac- 
companying pressure gradients which in turn drive large 
flow magnitudes (longitudinal, radial, elliptic) Q- 

The medium forms an "opaque core" which can stop 
energetic partons and which is surrounded by a less-dense 
"corona" from which single partons may escape radially 
outward while their in-going partners are stopped in the 
medium . That picture is said to account for sys- 
tematic modification of triggered dihadron correlations 
(usually with asymmetric pt cuts) [lO| and single-particle 
spectra (jet quenching and high-pt suppression) |ll| . 

However, recent analysis conflicts with the existence 
of a thermalized dense partonic medium and presents a 
different picture of A-A collisions. Untriqgered angular 
correlations (no pt cuts applied) [H, [H, [ij, [l^ [l^ re- 
veal that far from being "quenched" or lost to thermaliza- 
tion essentially all initially-scattered partons survive to 
produce correlated hadrons, albeit jets are signiflcantly 
modified in more-central A-A collisions. Jet correlations 
are dominated by minijets [l7| . Strong back-to-back 
(mini)jet correlations are still observed for all collision 
systems . Two-component analysis of pt spectra re- 



veals that the hard component is quantitatively described 
by pQCD for all A-A centralities, and the full parton 
spectrum is manifested as correlated fragments in the fi- 
nal state [H, [1^. The same hadron spectrum analysis 
also indicates that the scattered parton spectrum termi- 
nates near 3 GeV, not 1 GeV, implying a factor 25 less 
parton production (less density) than the conventional 
picture. 

Analysis of dihadron correlations typically incorpo- 
rates several important assumptions: 1) background sub- 
traction is determined in part by zero yield at minimum 
(ZYAM) dOl. 2) A contribution from "elliptic fiow" is 
determined by independent measurement of an azimuth 
quadrupole component in the form w| [2l|. 3) pt cuts 
define correlations from "hard scattering" as opposed to 
soft processes which could contaminate jet correlations. 
Certain features of dihadron correlations and interpreta- 
tions which support parton energy loss in a dense, ther- 
malized medium depend critically on those assumptions. 

It is important to resolve the striking incompatibility 
between two collision models emerging from the same 
underlying RHIC data. In this analysis I re-examine di- 
hadron correlation analysis methods, with special atten- 
tion to the ZYAM assumption, V2 subtraction and the 
consequences of trigger/associated pt cuts. 



II. DIHADRON AZIMUTH CORRELATIONS 

Modification of parton scattering and jet formation in 
more-central A-A collisions is seen as a possible signal of 
QGP formation. Jets are observed as per-trigger azimuth 
correlations by combinatoric reconstruction. Background 
subtraction, including ZYAM and V2 estimation, plays a 
central role in jet correlation analysis and interpretation. 



A. Dihadron correlation analysis 

High-pt dihadron azimuth correlations are used to 
study possible jet modifications in A-A relative to p-p 
collisions. The density of particle pairs which satisfy 
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certain pt (transverse momentum) cuts ("trigger" and 
"associated" particles) is plotted on azimuth difference 
A0 = ^associated " ^trigger- Trigger particlcs fall within 
one of two (usually asymmetric) pt intervals. Combina- 
toric jet reconstruction is a proxy for true event-wise jet 
reconstruction in A-A collisions with large multiplicities. 

A combinatoric background must be subtracted from 
the "sibling" (same-event) pair density to obtain the true 
distribution of jet-correlated pairs. The background is 
assumed to be modulated by an elliptic flow (W2) con- 
tribution. V2 is estimated by separate W2-specific anal- 
ysis. The zero reference is determined conventionally 
by the ZYAM assumption. In p-p and peripheral A-A 
collisions the dominant structures are a same-side (SS, 
|A0| < 7r/2) peak (nominally Gaussian), an away-side 
(AS, \A(j) — 7r| < tt/2) peak, and a cos(2A(^) sinusoid 
measured by w| (azimuth quadrupole). 

The ZYAM subtraction applied to more-central A-A 
collisions produces apparently interesting structures in 
the away-side component (nominally back-to-back jets) 
interpreted by some in terms of "Mach cones," possibly 
the result of ener getic p artons interacting with a dense 
QCD medium p^. [23l. [24| . However, the supporting argu- 
ments for ZYAM are questionable, and values of V2 used 
in the subtraction may be substantially overestimated 
leading to distortions of true jet correlations. 



B. Jet phenomenology 

An initial analysis of jet-related azimuth correlations 
in 200 GeV 0-10% central Au-Au collisions compared to 
p-p collisions imposed a high-pt cutoff (2-4 GeV/c) on the 
associated-particle interval as a means to insure that true 
jet phenomena dominated the correlations ^10^. Combi- 
natoric background subtraction invoked ZYAM and uti- 
lized values of V2{pt) from independent measurements. 
Mean value V2 — 0.07 was used for the subtraction. The 
correlation analysis followed observation of high-pt sup- 
pression in single-particle spectra [ill • The major conclu- 
sion was "disappearance of the away-side jet" (apparent 
back-to-back jet suppression). 

Followup studies 2^ [2^ sought to determine where 
the energy "lost" at larger pt appears in the final state. 
Jet correlations below 4 GeV/c were said to be af- 
fected by hydrodynamic flows and quark coalescence be- 
lieved to dominate hadron production in that pt inter- 
val [2^ . Combinatoric techniques permitted extension of 
associated-particle Pt cuts to as low as 0.15 GeV/c. 

Azimuth correlations from conventional (ZYAM) anal- 
ysis have the following characteristics: When smaller-pt 
associated particles are included both SS and AS peaks 
are broadened on azimuth in more-central Au-Au colli- 
sions, and the AS pt spectrum is "softened" [2^. The SS 
peak is enhanced relative to p-p, and the AS peak devel- 
ops a minimum (dip) at A(/) = tt. Those systematics are 
thought to reveal "characteristics of the energy transport 
of the quenched partons" on both pt and A(/) ^25^1 . 



C. Conventional physics interpretation 

The results of [l3| were interpreted to reveal an opaque 
"core" (dense medium) formed in more-central Au-Au 
collisions which stops most partons produced in or pass- 
ing through that region. The core was interpreted as a 
thcrmalizcd QCD medium with large matter and energy 
densities ^, i9(|. Partons which interact strongly in the 
medium lose most of their energy by radiating gluons 
and are thermalized (25j . "The away side jet traverses a 
large amount of matter" [emphasis added] |2fl] . "Shower" 
gluons may be emitted at large angles, possibly resulting 
in novel jet structures [23l . [2^ . pQCD calculations of 
parton energy loss require 30 times normal nuclear gluon 
density to account for central Au-Au results [27l |. 

In the core-corona model of A-A collisions some par- 
tons are produced in a less-dense "corona" region out- 
side the dense core and may survive to fragment outside 
the A-A system |25|]. High-p^ hadrons and hadron pairs 
mainly come from such partons, which suffer minimal in- 
teraction with the medium. Partons with pt > 5 GeV/c 
are said to exhibit in-vacuum fragmentation based on 
1) a p-p- like single-particle power-law spectrum and/or 
suppression (e.g., Raa [28]) independent of pt, and 2) 
p-p-like azimuth correlations with distinct SS and AS 
Gaussians. 

According to [2^ increased soft-particle production in 
Au-Au compared to p-p indicates parton thcrmalization. 
Final hadron remnants in central Au-Au collisions no 
longer exhibit jet-like correlations. Persistence of a broad 
AS peak structure at smaller pt indicates global momen- 
tum conservation. Interactions seem to drive particles 
from jet fragmentation and from the bulk medium to- 
ward mutual equilibrium, possibly implying a high de- 
gree of thcrmalization of the "medium." A significant 
amount of associated-particle energy may come from the 
medium through recombination, scattering or fiow. 



D. Alternative interpretations 

Recent analysis suggests alternative interpretations of 
spectra and angular correlations. Two-component anal- 
ysis of 200 GeV Au-Au spectra [H, [3 suggests that the 
spectrum hard component represents fragments from the 
entire initial-state parton spectrum. Fragmentation func- 
tions are modified in more-central A-A collisions, but no 
partons are lost to thcrmalization. Analysis of 2D angular 
correlations [Till is consistent with that picture: minijet 
correlation strength is proportional to N-N binary col- 
lisions (niiinary) for peripheral collisions, but increases 



more rapidly than nu^ 



for more-central collisions. 



Triggered azimuth correlations are not easily compared 
with normalized per-participant minimum-bias results. 
Such comparisons could serve to test the extent of par- 
ton thcrmalization by direct comparison to pQCD predic- 
tions 19]. High-pt trigger particles are themselves sub- 
ject to changes in A-A dynamics with centrality. In con- 
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trast, absolute spectrum and correlation measures reveal 
that while AS correlations are strongly altered, scattered 
partons are not thermalized (stopped in a medium). 

The SS jet width in ID azimuth correlations does not 
increase significantly with A-A centrality as would be 
expected from strong parton interactions with a dense 
medium (multiple scattering or gluon bremsstrahlung). 
In fact, for minimum-bias 2D correlations the SS peak 
azimuth width is observed to decrease with increasing 
centrality No SS jet structure corresponds to "Mach 
cones." Persistence of a broad AS peak at smaller pt has 
been attributed to global momentum conservation [26| . 
but the AS peak amplitude scales with binary collisions, 
not with participant pairs as would be expected for global 
conservation. 

Imposition of a 2-4 GeV/c associated cut in led 
to the conclusion that the AS jet is eliminated in more- 
central A-A collisions, implying absorption (thermaliza- 
tion) of one of the scattered partons. However, rejecting 
associated particles in the 0.15-2 GeV/c interval intro- 
duces a large bias. Extension of the AS pt interval down 
to 0.15 GeV/c reveals copious AS jet correlations. AS 
correlation structure is dominated by interjet or jet-jet 
(parton-parton) correlations, not intranet correlations. 
AS peak modifications per se may not reflect changes 
in parton abundance or fragmentation. 

Effects of kt broadening may dominate AS correlation 
structure. The role of initial-state parton multiple scat- 
tering in A-A collisions is still not fully appreciated, al- 
though kt broadening has been studied with Monte Carlo 
simulations (loj . The role of independently-measured V2 
in combinatoric background subtraction which isolates 
jet correlations is also not well-defined. Oversubtraction 
of V2 can (and does) result in substantial artifacts. 

Production of a dense QCD medium could indeed re- 
sult in jet quenching and other QGP "signals," but modi- 
fications to parton scattering and fragmentation observed 
to date don't require production of a QGP. Given the im- 
portance of angular correlations and jet-related phenom- 
ena to interpretation of RHIC data we should reconsider 
the conventional picture in light of recent results. 



III. THE ZYAM SUBTRACTION PROCEDURE 

The main object of "triggered" dihadron correlation 
analysis is systematic study of jet structure variations 
with A-A centrality and trigger/ associated combinations 
of Pt cuts. Accurate isolation of jet correlations from 
other correlation structure is technically challenging. 
The distribution of "sibling" pairs (pairs within single 
events) averaged over an event ensemble contains contri- 
butions from several sources: 1) uncorrelated combina- 
toric background, 2) back-to-back jet pairs, 3) isolated 
jets (partner outside the rj acceptance), 4) an azimuth 
quadrupole component conventionally attributed to el- 
liptic flow and 5) small contributions from other sources 
(e.g., quantum correlations, electron pairs from gamma 



conversions) . 

To isolate jet contributions 2) + 3) in the ZYAM pro- 
cedure a subtraction is performed. A nominally un- 
correlated mixed-pair reference distribution 1) approx- 
imates the combinatoric background structure (e.g., ac- 
ceptance/efficiency azimuth variations). But normaliza- 
tion of mixed pairs relative to sibling pairs depends in 
part on fluctuation sources (e.g., impact parameter fluc- 
tuations, large-scale correlations) which are not reliably 
estimated at present. The normalization is instead de- 
termined by the ZYAM criterion. The quadrupole com- 
ponent 4) is estimated from measurements of V2{pt) by a 
variety of methods which are susceptible to various "non- 
flow" errors (mainly jet correlations). Background 1) is 
then modulated by the estimated quadrupole amplitude, 
assuming a thermalized flowing bulk medium. 

A. Dihadron correlations 

Absent particle identiflcation two-(chaTged)-particle 
correlations are measured. Electron pairs due to pho- 
ton conversions from tt^ decays appear near the angular 
origin, are prominent in 2D angular autocorrelations [l3| 
and may distort the ID same-side (jet) peak on azimuth. 
The electron-pair structure is narrow on pseudorapidity 
as well as azimuth difference, and is thus suppressed to 
some extent in the ID azimuth projection. 

Dihadron correlations have three main components: 
1) a same-side (SS) peak centered at the origin, nom- 
inally Gaussian in shape and representing intrajet cor- 
relations, 2) an away-side (AS) peak centered at tt and 
representing internet correlations from back-to-back jets, 
3) a cos(2A(/)) term conventionally measured by The 
cos(2A(/)) term is the azimuth quadrupole component [29| . 

B. Trigger /associated pt cuts 

"Triggered" jet correlation analysis emulates event- 
wise jet reconstruction in the high-multiplicity environ- 
ment of A-A collisions by imposing pt cuts on particle 
pairs. A trigger particle at larger pt is intended to esti- 
mate the leading-parton momentum. An associated par- 
ticle at smaller pt should then sample the jet fragment 
distribution. Aside from that model-dependent language 
the resulting azimuth correlations correspond to speciflc 
cut boxes on the space {pti,Pt2)- Minimum-bias or "un- 
triggered" jet correlations are obtained by applying no pt 
cuts other than an acceptance cut near 0.1 GeV/c. 

An asymmetric cut system has some validity for trigger 
Pt > 4: GeV/c where the jet multiplicity is signiflcantly 
larger than 2. For the majority of jets (minijets) that con- 
dition is not satisfied. A symmetric pt cut then minimizes 
bias [isj . Asymmetric pt cuts imply significant distance 
from the diagonal of the (j)ti,Pt2) cut space, which has 
implications for joint V2(j>ti,Pt2) distributions compared 
to marginal distributions and factorization assumptions. 
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C. V2 estimation 

Azimuth correlation measure V2 is estimated by sev- 
eral methods denoted by W2{inethod}, e.g. V2{2} (two- 
particle correlations), V2{EP} (event plane or "stan- 
dard" method), W2{4} (four-particle cumulants) , etc. [soj . 
In practice, different methods return different V2 val- 
ues for the same physical conditions, depending on the 
amount of "nonflow" (minijets) confused by each method 
with the intended quadrupole component (29l . [sij . 

Various strategies are employed to reduce contribu- 
tions from nonflow and V2 fluctuations. V2{4} is said 
to eliminate nonflow contributions assuming nonflow 
is "clusters" of 2-3 particles (inconsistent with minijet 
systematics in more-central Au-Au collisions). Aver- 
age {v2{2} + w|{4})/2 is said to eliminate contributions 
from V2 fluctuations based on several physical assump- 
tions [3l|. In the simulations presented here w|{2} is 
determined by a fit (light dotted curve) to the entire az- 
imuth distribution, and the background quadrupole am- 
plitude V2 is defined by u|{2}/2 (assuming w|{4} = for 



central A-A collisions) 



is in that case locked to the 



amplitude of the SS Gaussian. 

The quadrupole contribution relevant to two-particle 
correlations subject to pt cuts is determined by joint 
distribution w|(pti:?'t2), which has not been measured. 
In the conventional approach factorization is assumed — 
V2iPti,Pt2) ~ V2{pti)v2{pt2), where V2{pt) is from a 
separately-measured marginal distribution. 



D. ZYAM background subtraction 
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FIG. 1: Left: Simulated "raw" (unsubtracted) dihadron 
correlations (solid dots), azimuth quadrupole amplitude A2 
with corresponding sinusoid (light dotted curve) and ZYAM 
subtracted background (bold dotted curve). The ZYAM- 
estimated background offset is the dotted line. Right: Re- 
sult of ZYAM background subtraction in the left panel. The 
minimum at n in the AS (away-side) peak is notable. 



since wfl^} = holds for R with no quadrupole compo- 
nent. The bold dotted curve is B from Eq. with Bq 
adjusted so that difference S = R — B satisfies ZYAM. 

Fig. [T] (right panel) shows resulting "jet" correlations 
S as the points. The "raw" jet component has been sub- 
stantially reduced in amplitude, and the AS peak has a 
minimum at tt. That is the general form of dihadron 
correlations for more-central A-A collisions after ZYAM 
subtraction. The solid curve is a free fit to the distri- 
bution discussed further below. The systematic error of 
the subtraction is sometimes estimated by the statistical 
error in a small interval about the minimum (e.g. psj]). 



The per-event and per-trigger "raw" pair density is 

R{A(b;pn,Pt2) = ^,5(0, - 0, - A(/.) \(1) 

averaged over an event ensemble and integrated over 
pseudorapidity acceptance A77, where pti and pt2 denote 
"trigger" and "associated" cut bins and Ntrig = X^iepti ' 
By hypothesis R includes jet-correlated pairs S (the jet 
signal) and (uncorrelated -I- V2) background B. To iso- 
late jet correlations S background B must be subtracted. 
The background form is assumed to be 

B{A^;pti,pt2) = Bo + B2ipti,Pt2) cos{2A^) (2) 
= Ba[l + 2v2ipti)v2iPt2) cos{2A(j))], 

where coefficient Bq is determined by ZYAM and 
'L>2{pti) 'L>2ipt2) V2 is the actual value estimated from 
one or more w|{mcthod} measurements. 

Fig. [T] (left panel) illustrates ZYAM subtraction. 
The points simulate a measured R distribution (with 
quadrupole component = 0). The light dotted curve rep- 
resents w|{2} obtained by fitting the entire distribution 
with cos(2A(/))-|-constant (consistent with the definition 
of u|{2}). The value of vl adopted for B is u|{2}/2. 



IV. AZIMUTH CORRELATION STRUCTURE 

In p-p collisions triggered ID azimuth correlations are 
approximated by Gaussian peaks centered at and tt ra- 
dians. In A-A collisions the structure evolves strongly 
with centrality. The AS peak broadens and a cos(2A(/>) 
(quadrupole) sinusoid dominates the structure for inter- 
mediate centralities. Either jet peak may deviate from 
a Gaussian shape. The main question is how to isolate 
true jet correlation structure while minimizing bias. 



A. Periodic peak arrays 

Distributions on azimuth are periodic. Peaks (nomi- 
nally Gaussian) at (SS, same-side) and tt (AS, away- 
side) are elements of periodic arrays. The SS array is 
centered on even multiples of tt, the AS array on odd 
multiples. Nearby elements of an array outside a 27r in- 
terval can have a significant effect on distributions within 
the interval and should be included in fit models. Fail- 
ure to include nearby members of both peak arrays may 
result in significant fitting errors, especially for broader 
peaks. 
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A peak array (SS or AS) can be represented by a 
Fourier series of the form 

S'(A(?!);crA0,n) = Ao,„ + Ai,„ {1 + cos(A0 - n 7r)}/2 

OO 

+ ^ A,„^„ cos(m [A0 - n7r]), (3) 

m=2 

where the Am,n are functions of Gaussian width cta^, n 
is even for SS peak arrays (+), and odd for AS arrays 
(— ). The terms represent 2m poles, e.g. dipole (m = 
1), quadrupole (to = 2), sextupole (to = 3). As peak 
width (Ta0 increases the number of significant terms in 
the series decreases. The hmiting case is cta^ ~ 7r/2, 
for which a peak array is approximated by a constant 
plus single dipole term — [1 + cos(A0)]/2 (SS) or [1 — 
cos(A(/))]/2 (AS). For Gaussian peak arrays with cta^ < 1 
(and therefore nonzero sextupole) a Gaussian function is 
the more efficient representation. 

Fig. [5] (left panel) illustrates the sum of peak arrays 
(solid points) for SS and AS peaks extending beyond one 
27r interval. The SS Gaussian peak array is the dash- 
dotted curve, the AS array with cta^ ~ 7r/2 is the dashed 
curve (approximately pure dipole). The dotted curve is 
the quadrupole term of the SS array, which would add a 
large "nonflow" contribution ^2,+ to vll-^} inferred from 
that distribution (see next subsection). 
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FIG. 2: Left: Periodic arrays of SS (same-side, dash-dotted) 
and AS (away-side, dashed) peaks. The SS peaks are Gaus- 
sians. The AS peaks are well-described by a dipole. The 
dotted sinusoid corresponds to the m — 2 Fourier component 
of the SS peaks. Right: Fourier amplitudes Am of a Gaussian 
[Eq. ([Sjl] vs peak width aA</>. 

Fig. [2] (right panel) shows the Fourier amplitudes of 
a Gaussian peak for the first five terms of Eq. ^ as 
functions of the r.m.s. peak width. For cta^ ^ n/2 only 
the constant and dipole terms contribute. For narrower 
peaks terms with to > 1 become significant. 

B. The quadrupole component 

Analysis whose main objective is "elliptic flow" mea- 
sures some part of the total quadrupole component of az- 
imuth correlations depending on the specific "method." 
There are three phenomenological sources of the total 



quadrupole component: the SS jet peak, the AS jet peak 
and the non-jet (NJ) background. The NJ component 
may or may not be related to a hydro phenomenon. 

The m — 2 Fourier amplitude of the entire azimuth 
distribution is A2 — 2Aowf{2}, based on the defini- 
tion of u|{2} ~ u|{£'P}. The to = 2 Fourier terms 
from the SS and AS peak arrays contribute to ^ll^} 
a combined "nonflow" systematic error 62- If the AS 
peak is broad (typical for more-central A-A collisions 
and/or untriggered correlations) "nonflow" 62 is dom- 
inated by the SS peak. The following correspondence 
can be made to conventional terminology and symbols. 
A2^± = 2ylo,±^2,±- ^2,-1- + ^2,- — ^2 (conventional 
"nonflow" measure). vl{EP} « w|{2} ^62+ v'^{NJ}. 
w|{2Z)} ~ 82 - -l-u|{iVJ} (from untriggered 2D autocor- 
relations). For typical untriggered correlations 82.- ^ 
(^2,+ , f|{A^J} and V2{2D} « V2{NJ}. Background com- 
ponent B2 — 2 Bq V2^NJ} may or may not be caused by 
a hydro mechanism [31( |. 



C. Peak models 

Concern has been expressed that the background off- 
set cannot be determined accurately by peak model fits 
if the exact form of the jet peaks is not known a pri- 
ori. It is claimed that due to in-medium modifications 
peak shapes may not be Gaussian. "Rigorous decom- 
position of the jet from its underlying event currently 
requires assumptions about the jet shape or the physics 
of the underlying event" [H]. Studies with PYTHIA are 
interpreted to conclude a substantial difference between 
peak shapes with and without radiation effects, including 
non-Gaussian shapes. Given such difficulties the ZYAM 
procedure is invoked as an alternative. 



■a 




FIG. 3: Simulated dihadron correlations for p-p collisions 
from PYTHIA [25| with cuts 3-5x3-5 GeV/c applied illustrat- 
ing deviations from Gaussian peak structures with (upper) 
and without (lower) initial- and final-state radiation invoked. 

Fig. m shows PYTHIA simulations (samples from Fig. 5 
of [2^) which reveal changes in peak shapes with and 
without initial- and final-state radiation. The simulation 
is invoked to support the conclusion that jet peak shapes 
are sufficiently uncertain (i.e., non-Gaussian) that peak 
models cannot be used to determine the background off- 
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set, that the alternative ZYAM method is necessary. 
To test that assertion the modificd-Gaussian form 

A 



1 



2n \^ 



(4) 



is compared to data. That expression is to Gaussians 
what the Levy distribution is to exponentials [s^ . The 
Levy distribution, the average of an ensemble of exponen- 
tial distributions on variable x with fluctuating slope pa- 
rameter X (e.g., mt,T for Maxwell-Boltzmann distribu- 
tions), has the form l/{l + x/X n)", where 1/n = aj^/X'^ 
represents the relative variance of X over the ensemble. 
Eq. (|4]) describes an ensemble of Gaussians with fluc- 
tuating widths crA0 and relative variance cr^ /ct^ — 



are exp( 



x/X) and 



l/4n. Limiting cases for 1/n 

exp{-(A0/(TA0)V2}. 

Figure H] shows the samples in Fig. [3] as points. The 
solid curves through points are Gss + Gas + constant. 
The G parameters are noted and background constants 
are indicated by the solid lines. The solid curves accu- 
rately represent the simulations. The background offsets 
are well-determined. The dashed curves omit the back- 
grounds. The dash-dotted curves include Gss and Gas 
with l/n — — pure Gaussians G with the same widths. 
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FIG. 4: Left: Data with no radiation from Fig. [3] (points) 
described by modified Gaussians as in Eq. Q plus a constant 
background (solid curve). The peak parameters are noted. 
Dashed curves are with no background. Dash-dotted curves 
are with no width fluctuations (l/4n — 0). Right: Same, but 
for data with radiation invoked in PYTHIA. 

Figure [5] shows the same elements on a linear format. 
Especially for the SS peaks the difference between true 
and modified Gaussian is subtle. 

The good agreement between Eq. and the PYTHIA 
simulations reveals that deviations from Gaussians are 
primarily due to width fluctuations. For the extreme case 
of Fig. [3] (high pt cuts, fiA^ ~ 0.15) deviations of the peak 
shape from a Gaussian in p-p collisions are obvious on 
a semilog plot. However, the peak shape is amenable to 
exact description with a simple extension of the Gaussian 
model function to include fluctuations. 

For cut systems typically encountered in dihadron 
studies at RHIC the SS peak width is more typically 
0.4 to 0.5 rather than 0.16 as in the PYTHIA study 
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FIG. 5: Same as Fig. Ubut plotted on a linear scale. 



Other things being equal that implies a ten-fold decrease 
in relative variance l/4n. SS peak shapes should then be 
indistinguishable from true Gaussians. In central Au-Au 
collisions fluctuations in the SS peak width should be re- 
duced flve-fold by central-limit averaging, since there are 
on average 30 ± 6 minijets in each collision [l^. The 
mean AS peak width increases substantially in more- 
central Au-Au collisions, further reducing the effective 
l/4n, but driving the peak shape toward a pure dipole 
shape. 

For p-p collisions and pt cuts below 6 GeV/c, espe- 
cially for minijet (minimum-bias jet) structure, and for 
more-central Au-Au collisions, the SS peak shape is accu- 
rately modeled by a Gaussian. The AS peak is accurately 
modeled by three terms of a Fourier series (or a Gaus- 
sian in some cases). Background subtraction with small 
systematic uncertainties can be achieved by peak fitting. 



D. 2D angular autocorrelations 

2D angular autocorrelations on difference variables 
Va = Vi ~ V2 and 0a = (f'l ^ 4'2 provide substan- 
tial additional information beyond the ID azimuth pro- 
jection. The typical correlation structure includes a 
2D SS peak, an AS peak on azimuth (AS ridge) inde- 
pendent of pseudorapidity difference, and an azimuth 
quadrupole [29l.[3]|. 

The SS 2D peak is unambiguously isolated by model 
fits because of its variation on pseudorapidity difference. 
The AS ridge and azimuth quadrupole are uniform on 
77a and thus project onto ID azimuth without loss of 
information. The AS peak is typically dominated by 
the dipole term. For more-central A-A collisions, and 
generally for untriggered correlations, the AS peak may 
be pure dipole, in which case it is orthogonal to the 
quadrupole component. For peripheral A-A or p-p colli- 
sions and higher-pt cuts an AS Gaussian is more efficient. 
For intermediate cases (cr^^ ~ 1.0 - 1.2) the AS dipole 
model may allow a small contribution 62.^ to measured 
quadrupole amplitude w|{2il'} which then establishes an 
upper limit on any independent (nonjet) quadrupole com- 
ponent w|{iVJ}. 

Fig. [S] (left panel) illustrates untriggered (no special pt 
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cuts) 2D angular correlations for 200 GeV mid-central 



Au-Au collisions The SS peak is a 2D Gaussian 

which can be reliably isolated from the other correlation 
structure. The remaining structure then consists of [1 — 
cos(0a)]/2 dipole and cos(20a) quadrupole components. 





FIG. 6: Left: Simulated 2D angular autocorrelation showing 
an r;-elongated SS peak, an AS ridge and quadrupole. Right; 
The 2D autocorrelation with SS peak subtracted and pro- 
jected onto (p^ (points). The result is a dipole (solid curve) 
plus quadrupole (dotted curve). 

Fig. [H] (right panel) shows the distribution in the left 
panel with the 2D Gaussian fit removed and the dif- 
ference projected onto azimuth (points). The AS peak 
width for untriggered correlations is typically a^^ ^ 7r/2, 
and the peak is essentially a dipole. The AS peak and 
quadrupole are therefore orthogonal, implying negligible 
covariances among the three model components: SS 2D 
Gaussian, AS dipole and quadrupole. 

2D autocorrelations can also be constructed for specific 
(symmetric or asymmetric) pt cuts, providing combina- 
torial triggered jet reconstruction with reduced ambigu- 
ity. 2D histograms provide the most direct and accurate 
access to joint distribution V2{pti,Pt2)- There is no ev- 
idence in 2D autocorrelations for the AS double-peaked 
structures resulting from ZYAM subtraction. 



V. ZYAM CRITIQUE 

The purpose of this analysis is to examine the validity 
of the ZYAM procedure and consequent novel structures 
in reported jet correlations. In this section we review the 
logic and consequences of ZYAM. 



A. Assumptions supporting ZYAM 

The background level (zero offset) of sibling-pair dis- 
tribution R in Eq. ([1]) is well-defined when SS and AS 
peaks are completely separated, as in p-p collisions with 
"high-pt" cuts. For other collision systems or lower pt 
cuts the peaks may overlap, and the background level is 
not then visually obvious. For overlapping peaks and at 
higher pt, especially in more-central A- A collisions, u| 
subtraction is also ambiguous. 



The ZYAM procedure explicitly assumes that sub- 
tracted background B emerges from a thermalized bulk 
medium. Almost all particles, otherwise uncorrelated, 
are then correlated with the reaction plane and described 
by motivating the structure of B in Eq. ([2]). The 
problem remains to determine the true background offset 
for azimuth correlations, the constant B^ in Eq. In 
principle, the background for sibling pairs accepted in cut 
box {pti,Pt2) should be the mean pair ratio (ni n-i) jfii n2 
times the mixed-pair background normalized to pair 
number fii n2, in which case Bq — {rii 712). 

However, covariance (ni 712)— ni n2 (a fluctuation mea- 
sure) contains contributions other than jet correlations. 
Positive contributions include centrality (6) fluctuations 
and correlations with characteristic lengths exceeding the 
angular acceptance. Negative contributions include pos- 
sible local anticorrelations and, most importantly, canon- 
ical suppression of fluctuations if a centrality condition 
is imposed on rich within the pair angular acceptance or 
some part of it. Thus, Bq cannot be determined directly 
from present fluctuation [{rii 712)] measurements. 

Alternatively, Bq could be determined directly from 
correlation structure by model (e.g., Gaussian) fits to the 
SS and AS peaks. The main argument against model fit- 
ting is that the collision process could modify peak struc- 
tures away from a Gaussian shape. Incorrect fit models 
might then lead to biased background subtraction. If SS 
and AS peaks are overlapping ZYAM is invoked by that 
argument to determine Bq. 

ZYAM implicitly assumes that real distributions do 
not contain overlapping peaks, which contradicts many 
examples in nature including nuclear collisions. For most 
of the centrality range of Au-Au collisions and some 
Pt cuts the same-side and away-side peaks do overlap 
strongly, and the true background offset can differ greatly 
from the ZYAM estimate. 



B. Estimating vi(pti,pt2) 

The subtracted background includes w| correlations 
(azimuth quadrupole). In Eq. ^ joint distribution 
V2{Pti,Pt2) has been replaced by product ^2(^*1) ^^2(^*2) 
assuming factorization. The factors (marginal distribu- 
tions) are determined by separate analysis. Factoriza- 
tion is defended by the claim that in a thermalized bulk 
medium particles have no correlations among themselves, 
only with the common reaction plane. 

V2{pt) is measured by several methods denoted by 
W2 {method}, e.g., W2{2} (two-particle correlations), 
V2{EP\ (event plane), U2{4} (four-particle cumulants), 
etc. [30]- At larger W2{2} or V2{EP} is favored because 
of smaller particle yields. Both are strongly sensitive to 
minijets (nonflow). Extraction of an azimuth quadrupole 
component from 2D angular correlations provides alter- 
native measure W2{2-D} which is by construction inde- 
pendent of jet correlations to good approximation [31]. 

In conventional flow analysis it has been assumed that 
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w|{2} = w|+tT^2 and i;|{4} = v\—(j1^ , so that mean value 
(w|{2} + w|{4}) /2 = V2 eliminates "flow fluctuations" 
0' • The mean value is often used to determine 
V2{pt) factors for Eq. But v'i{2} = v^{2D}+62 which 
defines "nonflow" S2, dominated by the m = 2 Fourier 
component of the SS Gaussian (mini)jet peak (Blj . And 
f2{4} may also include minijet contribution S2 in more- 
central Au-Au collisions, because minijets then have mul- 
tiplicities > 4. Thus, (i;|{2} + v|{4})/2 ^ vl{2D} - 
+ (^2 + '52)/2, and independent (nonjet) quadrupole 
component d|{2_D} may be a small fraction of the non- 
flow (jet) terms for more-central Au-Au collisions. 

In Fig. [1] (left panel) the light dotted curve represents 
ti2{2} obtained (per its definition) by fitting the entire 
azimuth distribution with cos(2A0). By construction 
the data contain no independent quadrupole component. 
Thus, f2{2} — 82 is entirely "nonflow" or jet correlations. 
In essence, the quadrupole component of jet correlations 
is subtracted from jet correlations to produce distortions. 

C. A typical ZYAM analysis 

Figure [7] (left panel) shows simulated azimuth corre- 
lation data (points) for central [h = 0) Au-Au colli- 
sions. The distribution consists of SS Gaussian and AS 
dipole only. A fit to unsubtractcd data using a Gaus- 
sian-f dipoleH-quadrupole-fconstant model (solid curve) 
by construction returns the simulation parameters. Also 
shown is a fit to the distribution with Aq -I- A2 cos(2A0) 
only (dotted curve), which for this case corresponds to 
A2 = 2Aq vl{2} = 0.12 ~ P2/4 by definition of z;|{2}. 




2 4 2 4 

A(|) A^) 

FIG. 7: Left: Simulated azimutfi correlation data (points) 
including a SS Gaussian (dash-dotted curve) and AS dipole 
(dashed curve). The bold solid curve is a free fit to the data 
which returns the model parameters: ofi^set (PI), SS Gaus- 
sian amplitude and width (P2, P3), dipole amplitude (P4), 
quadrupole amplitude (P5). Right: The same data fitted with 
offset and SS Gaussian (as before), AS Gaussian amplitude 
and width (P4, P5) and quadrupole (P6). 

Figure [7] (right panel) shows a fit to the same "data" 
with two Gaussians and a quadrupole. In comparing the 
two panels several aspects are notable. 1) The offset lev- 
els are the same. 2) The properties of the SS Gaussian 
agree to a few percent. 3) The amplitude of the AS peak 



agrees to a percent. 4) The inferred quadrupole ampli- 
tude (zero in the simulation) is, in either case, 1% or less 
of the dipole amplitude. The inferred r.m.s. width 1.3 of 
the AS Gaussian is consistent with the dipole model. 

Figure [5] (left panel) uses the same simulated data to 
illustrate the ZYAM procedure. For this example ^ll^} 
is determined as above, and (ideal) {4} is zero. If the 
subtraction is defined by quadratic mean u| = (w|{2} -|- 
w|{4})/2 then since Ao ^ 8 the subtracted value is z;^ = 
0.12/(2 X 2 X 8), or V2 = 0.06 (compared to 0.07 in n^). 




2 4 2 4 



FIG. 8: Left: Simulated data as in Fig. [7] with ZYAM back- 
ground model (bold dotted curve). Right: Result of subtract- 
ing ZYAM background model from simulated data (open sym- 
bols). The bold solid curve is a free fit to the subtracted data 
which returns the original SS Gaussian and AS dipole model 
parameters plus the ZYAM offset and quadrupole parameters. 

Figure [8] (right panel) shows the result of ZYAM sub- 
traction as the points. The dash-dotted and dashed 
curves show the simulation input, the "right answer." 
The bold solid curve through the points is a free fit 
with the SS Gaussian-|-AS dipole-l-quadrupole model. 
The original parameters are returned in addition to the 
ZYAM-imposed w| and offset. In other words, no dis- 
tribution information is lost by the ZYAM subtraction, 
but the result is visually misleading. In particular, the w| 
oversubtraction results in a minimum in the AS peak, and 
the ZYAM offset subtraction imposes a large reduction 
in apparent jet yields. The same fit model was applied 
in Fig. [1] (right panel), with similar results. 



D. Consequences of ZYAM 

ZYAM subtraction with overestimated V2 values re- 
sults in distorted azimuth correlations, including a mis- 
leading double-peaked away-side structure convention- 
ally interpreted in terms of "Mach cones" [H, [l^l . Over- 
subtraction of V2 occurs for three reasons: 1) z^ilZ} 
(or u|{i?P}) includes large "nonfiow" (minijet) contri- 
butions for all A-A centralities but especially for more- 
central collisions, 2) w|{4} can also include minijet con- 
tributions for more-central collisions and 3) factorization 
of the joint w|(pti,j5t2) distribution can be questioned 
given asymmetric (off-diagonal) pt cuts. 

Fig. [5] (right panel) shows nominal jet signal S — R — 
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B as the points, illustrating the general form of more- 
central dihadron correlations after ZYAM subtraction. 
The original jet correlations are strongly reduced, and 
the away-side peak has a minimum at tt. The apparent 
reduction of jet correlations appears to confirm strong 
jet quenching and parton thermalization. The AS peak 
structure also suggests parton interaction with a dense, 
thermalized medium leading to shock-wave formation. 

The form of the ZYAM-subtracted "data" in Fig. [8] 
(right panel) depends only weakly on collision conditions, 
because w^l^} is dominated by the m = 2 Fourier com- 
ponent of the SS peak. Its magnitude is therefore locked 
to the jet-peak amplitude, and all correlation structure 
scales up and down together with collision centrality. 

Figure [S] (left panel) demonstrates that maxima near 
A(j) « TT ± 1 resulting from ZYAM subtraction are an 
inevitable consequence of V2 oversubtraction for typical 
angular correlations from RHIC collisions. For a generic 
SS/AS peak combination in more-central A-A collisions 
several values of A2 (integer multiples of 2Ao w| = 0.015) 
including zero (solid curve) are invoked in the subtrac- 
tion. The persistence of apparently-displaced peaks at 
the same locations is evident. Note the effect of ZYAM 
subtraction even for w| = (solid curve and dotted line) . 




A(|) A4) 



FIG. 9: Left: Demonstration of the effects of 112 oversubtrac- 
tion. Model jet correlations from Fig. [7] are the "data" (solid 
curve) from which quadrupole components of successively 
larger amplitudes are subtracted (dashed curves). Dotted 
lines represent zero levels established according to the ZYAM 
principle. Vertical solid reference lines are at tt ± 1.1. Right: 
ZYAM-subtracted distribution from Fig. |5] (right panel) fitted 
with single SS Gaussian and two AS Gaussians located sym- 
metrically about TT. Common AS Gaussian amplitudes and 
widths are P4 and P5. Displacements from n are P6. 

Figure [9] (right panel) illustrates attempts to charac- 
terize the AS peak structure with two Gaussians [1^. 
Nominal medium-induced "shoulder" Gaussians are lo- 
cated near A(j) w tt ± 1. In some cases a fragmentation- 
related "head" Gaussian is introduced at A(f> = tt. The 
bold solid curve is a free fit with shoulder Gaussians al- 
lowed to vary in centroid, width and amplitude, but sym- 
metrically about A(f> = TT. The detailed model fit can be 
contrasted with the light (not bold) solid curve of Fig. [8] 
(right panel), the "right answer" for these "data." 

In Figure [5] (right panel) the offset is allowed to vary, 
and the SS Gaussian amplitude is 0.36, closer to in- 



put value 0.50. If the offset were fixed at zero the SS 
peak amplitude would be 0.29, as in Fig. [S] (right panel). 
If the correct fit model were used the original simula- 
tion parameters would be returned. This exercise reveals 
the typical magnitude of systematic errors resulting from 
ZYAM and u| oversubtraction. 



E. Systematic errors and uncertainties 

ZYAM subtraction uncertainties derive from uncer- 
tainties in offset Bq and t;|. The ZYAM offset uncer- 
tainty is assumed to result only from determining the 
effective minimum of the unsubtracted (raw) pair distri- 
bution, which is in turn assumed to be dominate d by the 
statistical error within a small azimuth interval [25l. [2^. 

The present analysis demonstrates that the ZYAM as- 
sumption is contradicted by any distribution with over- 
lapping peaks. The true ZYAM systematic uncertainty 
is indicated by the typical difference between a ZYAM 
offset value and that inferred from peak fitting as in this 
analysis. The error can be large (up to 100% of true jet 
correlation amplitudes), leading to substantial underes- 
timation of jet correlations and related hadron yields. 

w| subtraction uncertainties are dominated by mea- 
surement uncertainties for marginal w|(pt) distributions. 
Especially for more-peripheral or -central A-A collisions 
minijet contributions to typically persist as large un- 
corrected systematic errors, i'l(pt) is estimated by aver- 
age (wi{2} + w|{4})/2 [11]. Where t;|{4} data are not 
available (e.g., 0-5% central Au-Au) the approximation 
w|{4} ~ t)|{2}/4 is sometimes made [2^. The minimum 
value of V2 (e.g., for central Au-Au collisions) is then 
'y|{2}/2 if (ideal) V2{4:} is known (i.e., zero) or 5'i;|{2}/8 
if f2{4} is not known. 

In central (6 = 0) A-A collisions w| = and ffl^} = 
62- Thus, the minimum value of w| assumed for back- 
ground subtraction in central collisions is 1 /2 the m = 2 
Fourier component of the SS jet peak. Jet correlations 
are converted to "nonflow" which is subtracted from jet 
correlations. Consequently, jet yields can be greatly un- 
derestimated. 



VI. ZYAM EXAMPLES FROM RHIC DATA 

Two examples from ZYAM-based jet correlation anal- 
ysis of RHIC 200 GeV Au-Au data are compared. ZYAM 
background subtraction is reversed to recover nominally- 
undistorted jet correlations. Free fits to data comparing 
a Gaussian-dipole-quadrupole model with a Gaussian- 
Gaussian-quadrupole model reveal consistent jet corre- 
lations in the two cases. 
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A. Dihadron correlations from STAR 

Fig. [TUl fupper left) shows data from Fig. 1 (upper left) 
of [11]. ZYAM-subtracted data from 200 GeV 0-12% cen- 
tral Au-Au collisions for (trigger x associated) 4-6x0.15-4 
GeV/c pt cuts are shown as solid points. Corresponding 
p-p data are shown as open circles. The bold solid curve 
is a free fit with offset (PI), SS Gaussian (amplitude P2, 
width P3), AS dipole (P4) and quadrupole (P5). The 
resulting offset is the solid line at Pi, the SS Gaussian 
is the dash-dotted curve, the dipole is the dashed curve, 
and the (negative) quadrupole is the dotted curve. 

Fig. [TO] (upper right) shows a reconstruction of the orig- 
inal ( "raw" ) data distribution prior to ZYAM subtraction 
based on Pi and P5 . The relation of p-p to central Au-Au 
data is quite different. Both SS and AS peaks increase by 
a factor six from p-p to central Au-Au, a large increase 
in the jet yield which is not apparent in the upper-left 
panel determined by conventional ZYAM subtraction. 
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FIG. 10: Upper left: ZYAM-subtracted angular correlations 
for 0-12% central 200 GeV Au-Au collisions (solid points) ^ 
with free fit (bold solid curve) of SS Gaussian (dash-dotted 
curve), AS dipole (dashed curve) and quadrupole (dotted 
curve). Open points are p-p data relative to ZYAM zero. 
Upper right: ZYAM subtraction reversed, true zero level re- 
covered from free fits to data (solid points), compared to p-p 
data (open symbols). Lower left: Consequences of raising the 
associated pt cut to 2 GeV/c [2^. Most of the SS fragment 
distribution is rejected and AS correlations are strongly sup- 
pressed. Lower right: Fit of two Gaussians plus offset to p-p 
data. 

Without the unsubtracted data we cannot determine 
the original quadrupole component exactly. However, 
from information in the upper-right panel we can esti- 
mate what value of 2 Aq u| was used in the subtraction. 



The SS peak amplitude is P2 = 3.37 and the width is 
P3 = 0.50. Fig. [2] (right panel) indicates that m = 2 
amplitude A2 for that width is 25% of the SS amplitude. 
Using the prescription = ('u|{2}-Hw|{4})/2 v^{2}/2 
for central collisions the minimum value of 2 Ao u| = 
0.25 X 3.37/2 = 0.42 from the definition of w|{2}. The 
corresponding number from ZYAM subtraction in the 
upper-left panel is —2 P5 = 0.42, suggesting that the 
quadrupole component of the unsubtracted data was con- 
sistent with zero. 

Fig . [TOl (lower left) shows data from Fig. 1 (lower left) 
of [26| with the same fit model applied but with the as- 
sociated Pt cut raised to 2-4 GeV/c. While the AS peak 
is strongly reduced by the elevated associated-pf cut it is 
still described by a dipole to the error limits of the data. 

Fig. [TO] (lower right) shows a two-Gaussian fit to the 
p-p data from [2^. Comparing p-p and Au-Au data the 
SS peak width changes only slightly from p-p (0.4) to 
central Au-Au (0.5). The AS peak width doubles from p- 
p (0.7) to Au-Au (~ 1.4). The latter should be expected 
due to increased fct broadening of the AS peak, since 
the mean participant-nucleon path length (as N-N binary 
collisions) increases from 1 to 6 in central Au-Au [34| . 

Figure 2 (upper panel) of [1^ reviews the apparent 
centrality systematics of SS and AS peak integrals. In 
either case an increasing trend for peripheral collisions 
appears to saturate for more-central collisions. However, 
from Fig. [To] (upper-right panel) we find that the 0-12% 
Au-Au/p-p ratio for both SS and AS is about six times 
that for p-p or peripheral Au-Au collisions — greatly ex- 
ceeding what is implied by the ZYAM subtraction. A 
similar reconsideration of [10| with reduced pt cuts might 
reveal that the "away-side jet" is still present, albeit with 
altered fragment (associated-particle) pt distribution. 



B. Dihadron correlations from PHENIX 

Figure [TT] (left panel) shows ZYAM-subtracted az- 
imuth correlation data from Fig. 9 of [2^ (points) corre- 
sponding to 2-3x2-3 GeV/c pt cuts. The data are from 
0-5% central Au-Au collisions at ^/sWn = 200 GeV. The 
bold solid curve is a free fit with offset (PI), SS Gaussian 
(amplitude P2, width P3) and away-side Fourier series 
(dipole P4, quadrupole P5, sextupole P6). The dash- 
dotted curve is the fitted SS Gaussian, the dashed curve 
is the fitted dipole, and the dotted curve is the fitted 
negative quadrupole. 

Figure [TT] (right panel) shows a fit replacing the AS 
Fourier series with an AS Gaussian (amplitude P4, width 
P5) plus quadrupole (P6), which describes the data 
equally well. The SS peak parameters are the same in the 
two panels. In the left panel the Fourier series describing 
the AS peak is dominated by dipole and quadrupole. The 
sextupole Pg, which would be significant if an AS Gaus- 
sian were required, is negligible. The quadrupole is large 
and negative, reflecting mainly what was subtracted in 
the ZYAM procedure. The free fit recovers the correct AS 
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FIG. 11: Left: ZYAM-subtracted correlation data (open 
circles) with free fit (bold solid curve) of SS Gaussian 
(dash-dotted curve), AS dipole (dashed curve) and negative 
quadrupole (dotted curve) plus offset (solid line). Right: The 
same data fitted with SS Gaussian (P2, P3), AS Gaussian 
(P4, P5) and quadrupole (P6), leading to equivalent results. 
The solid points are the open points translated by offset Pi. 



jet structure and background offset. In the right panel 
the AS Gaussian contains the same information as the 
Fourier terms in the left panel. Gaussian width 1.3 is 
consistent with a pure dipole representation. 

The quadrupole component of the unsubtracted data 
can be estimated for this case as well. The SS peak 
amplitude is P2 = 0.056 and the width is P3 = 0.55. 
The 2Aov^ estimate is then 0.25 x 0.056/2 = 0.007. 
The quadrupole amplitude after ZYAM subtraction is 
— 2P5 = 0.01, consistent with w| 3w|{2}/4 and negli- 
gible net quadrupole component in the unsubtracted data 
(likely for 0-5% central data). 

Fig. [TT] reveals that the increase in the SS peak yield 
over ZYAM subtraction is (0.056x0. 55)/(0.024x0.3) - 4, 
and the increase in the AS yield is (0.038 x 1.3)/(2 x 
0.008 X 0.3) ^ 10. The increase can be appreciated vi- 
sually in the right panel by comparing the solid points 
to the bold solid curve. Again, jet yields are strongly 
suppressed by conventional ZYAM and V2 subtraction. 



VII. DISCUSSION 

I first discuss problematic aspects of conventional anal- 
ysis and interpretation of triggered azimuth correlations 
prior to background subtraction. I then consider the con- 
sequences of ZYAM and v2 oversubtraction: underesti- 
mation of jet yields and distortion of jet correlations. 
Finally, I discuss contrasts with untriggered angular cor- 
relations and a proper fit model. 



A. The myth of the "away-side" jet 

In the context of so-called triggered jet analysis con- 
ventional language refers to an "away-side" (AS) jet [l^, 
[25l . [2^ . AS azimuth correlation structure is attributed to 
single jets which are complementary to "high-pt" trigger 



particles (i.e., proxies for trigger partons). The AS jet is 
said to broaden or, in the case of double-peaked struc- 
ture, to manifest gluon bremsstrahlung radiation or Mach 
shocks [23]. The description implies that two classes of 
jets exist: SS jets from the "corona" which are relatively 
intact (in- vacuum), and AS jets which suffer strong inter- 
actions with a dense medium (core) , supporting the core- 
corona collision model. That language seems to represent 
a misunderstanding of jets in A-A collisions. 

AS azimuth structure represents internet correlations, 
including not only znirajet structure (individual parton 
fragmentation) but parton-parton angular correlations as 
well, especially kt contributions (e.g., dijet acoplanarity) . 
kt generation by initial-state parton multiple scattering 
in A-A collisions may dominate AS correlations there. 

The SS jet peak on azimuth reveals mfrajet structure 
and represents all jets within an acceptance. Especially 
for untriggered angular correlations any modifications to 
jet structure (including Mach cones) should be mani- 
fested in the SS peak. Changes in AS correlations rel- 
ative to the SS peak should first be attributed to the 
parton-parton relationship, dominated by kt production 
during initial-state N-N (parton-parton) scattering. 

The role of kt (nucleon-intrinsic and from initial-state 
parton multiple scattering) even in p-p collisions is quite 
significant. The AS peak width in p-p collisions is typi- 
cally 50% larger than the SS peak, reflecting jet acopla- 
narity from the kt contribution to parton-parton trans- 
verse center-of-mass motion. The mean nucleon path 
length (as N-N encounters) is six in central Au-Au, im- 
plying up to six-fold increase in kt variance or 2.5-fold 
increase in AS peak width, which then accounts for typi- 
cal AS widths tTc6 ~ 1.3 — 1.5 in central Au-Au collisions. 



B. Consequences of trigger /associated pt cuts 

Application of pt cuts to emulate jet reconstruction 
can strongly bias jet-related azimuth correlations, leading 
to misinterpretations. To interpret the consequences of 
Pt cuts correctly a valid model of parton fragmentation 
(fragmentation functions or FFs in vacuum and medium) 
and parton-parton correlations {kt effects) is required. 

For single-parton fragmentation (e.g., "in- vacuum" 
jets) in 200 GeV p-p collisions the following trends are 
observed [H [fl [l|]: 1) The peak (mode) of FFs oc- 
curs in 1-2 GeV/c for all parton energies from 3 GeV 
up to several hundred GeV [s^. 2) p-p FFs are strongly 
suppressed below 1 GeV/c compared to jets in e'^-e~ col- 
lisions. 3) About half the fragments from minimum-bias 
jets or "minijets" in p-p collisions, with most-probable 
parton energy = 3 GeV, appear below 1 GeV/c. 

Some implications of those systematics are: 1) A typ- 
ical trigger pt cut at 4 GeV/c eliminates most scattered 
partons, since pt,parton > Pt,trigger, the mode of the par- 
ton spectrum is near 3 GeV and the spectrum is falling 
(xp^"^'^. 2) FF suppression at smaller fragment momen- 
tum in p-p collisions already narrows jet angular correla- 
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tions relative to in- vacuum jets from e-e collisions. 3) An 
associated cut at 2 GeV/c (e.g., [l3|) eliminates most 
of the FF for any parton energy, further narrowing jet 
angular correlations (only the "tip" of the jet survives). 

A trigger pt cut biases the underlying parton spectrum 
to larger energies and so could in principle be useful for 
systematic energy-loss studies. However, the trigger cut 
also biases the associated kt distribution, especially in 
more-central Au-Au collisions where accumulated kt can 
be large. Aside from jet (parton) acoplanarity {kt broad- 
ening of the AS azimuth peak) a unique consequence of 
the trigger cut is a kt bias in which the parton-pair kt sum 
is preferentially oriented nearly parallel to the trigger- 
parton momentum, thereby increasing the trigger-parton 
momentum and decreasing the partner momentum in the 
laboratory frame. Such kt bias can give an impression of 
large AS energy loss in a dense medium when the SS/AS 
momentum asymmetry is actually due to initial-state kt- 

In-medium modification of FFs in more-central Au- 
Au collisions can be described by alteration of the split- 
ting process [36j . wherein energy is not lost from a jet, 
is instead redistributed on pt to a degree depending on 
Au-Au centrality. More jet-correlated fragments appear 
at smaller pt [l9|. The effect is observed in untriggered 
two-particle correlations, including a sharp transition at 
intermediate Au-Au centrality between unmodified p-p 
jets and strongly- modified but still intact jets in more- 
central Au-Au [3, [3]. Modification of FFs in central 
Au-Au collisions can transport much of the jet below an 
associated-particle pt cut interval, giving the false im- 
pression that back-to-back parton pairs are suppressed, 
i.e., that one parton is "thermalized" in a dense medium. 



C. Consequences of ZYAM and V2 oversubtraction 

This analysis demonstrates that for any case where 
there is significant SS and AS peak overlap the ZYAM 
assumption discards a significant fraction of the true jet 
yield. In more-central collisions up to 90% of the jet yield 
may be discarded as a result of ZYAM plus u| oversub- 
traction, as in Fig. [11] (AS peak). Oversubtraction of u|, 
which is guaranteed by conventional subtraction meth- 
ods, further reduces the apparent jet yield and introduces 
strong distortions in AS correlations. The resulting dis- 
tortions can be misinterpreted as evidence that a large 
fraction of scattered partons is thermalized in a dense 
medium, and that novel QCD phenomena result. 

In [l^ the double-peaked AS structure induced in 
two-particle correlations is duplicated in a more-complex 
three-particle correlation analysis. The analysis is in- 
tended to demonstrate that relative to a trigger particle 
some pairs of AS particles are widely separated by a char- 
acteristic angle which might imply conical emission from 
a parton passing through a dense medium, as in Mach- 
cone shock waves [2^ [IJ . 

It was found that the distance of AS "Mach peaks" 
from TT ( "cone opening angle" ) does not depend on A- 



A centrality or associated-particle pt- The fixed open- 
ing angle has been interpreted to measure the speed of 
sound in the QCD medium [H, [IJ] . Yet the properties 
of a true QCD medium ought to depend strongly on A-A 
centrality (e.g., compared to p-p collisions). The present 
analysis indicates that the AS double peak is an artifact 
having a fixed form for simple algebraic reasons, with no 
relation to collision dynamics or a QCD medium. 

D. Untriggered 2D angular autocorrelations 

Untriggered 2D angular autocorrelations [H, [H, [3, 
[l5l . [l6j differ from conventional triggered dihadron 
correlation analysis in several ways: 1) Absent trig- 
ger/associated Pt cuts the entire fragment yield from the 
entire minimum-bias parton spectrum is measured. 2) 
No ad hoc (ZYAM) offset is imposed; the offset is in- 
ferred from free fits to correlation structure. 3) Model 
fits impose only minimal a priori assumptions which are 
subsequently tested with data. 4) No external V2 mea- 
surement is imposed; an independent quadrupole com- 
ponent is inferred directly from the model fit. 5) The SS 
peak is uniquely isolated in the 2D fit due to its rj depen- 
dence. 6) Remaining structure is described by model- 
independent azimuth Fourier amplitudes. The measured 
quadrupole amplitude serves as an upper limit to any 
nonjct quadrupole component. 

Unsubtracted dihadron correlations for pairs within 
specific Pt cut bins serve as the best measurement of 
joint distribution v\{pti,pt2) simultaneously with corre- 
sponding jet structure. ZYAM subtraction disregards 
direct i;| measurements, instead invokes products of 
V2{pt) marginal values of questionable accuracy (large jet 
or "nonflow" contribution likely). Conventional V2{pt) 
measurements for 0-5% central Au-Au collisions imply 
\/v2{pti) V2{pt2) ^ 0.06 — 0.07, whereas untriggered 2D 
angular correlations return V2 < 0.01. The difference in 
quadrupole amplitude ?;| is a factor 30-50. 

The SS 2D Gaussian is tightly constrained by the 2D 
fit. Its amplitude is not available, as in a ID fit, to offset 
distortions of the AS peak imposed by w| (quadrupole) 
oversubtraction. No double-peaked AS structure has been 
observed in 2D angular autocorrelations. 

E. Recommended procedures 

Whatever the subsequent analysis, unmodified di- 
hadron pair distributions with absolute pair numbers per 
collision event should be published for all collision sys- 
tems. A neutral presentation method for untriggered 
angular correlations normalizes the number of mixed 
pairs (reference) to the number of sibling (same-event) 
pairs [H, [l3|- The integral of the net-pair distribu- 
tion (number of "correlated" pairs) is made zero by con- 
struction to avoid suggesting that the net-pair number 
is known a priori. The true offset can be estimated sub- 
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sequently via model fits. Mean event multiplicities for 
accepted particles should be included in the presentation. 

The fit model should be as simple as possible (min- 
imum parameter number) yet physically interpretable. 
Peaks can be distinguished from background by model 
fits as long as the peaks are resolved according to the 
Rayleigh criterion (sum of r.m.s. widths significantly less 
than peak separation tt, which is typically satisfied for 
azimuth correlations). The Gaussian shape for the SS 
peak may be questioned in isolated cases (cf. Sec. lIVCI) . 
but for most cases of relevance to RHIC data the ID (on 
azimuth) or 2D Gaussian model should be adequate. 

For higher pt cuts and/or p-p collisions the AS azimuth 
peak may be described by a Gaussian. For untriggered 
correlations and/or for more-central Au-Au collisions the 
AS peak width increases substantially, the higher Fourier 
components of the AS Gaussian decrease, and only a few 
terms of a Fourier series may be required. In the asymp- 
totic limit (e.g., tr^ > 1.3) periodic AS peak structure 
simplifies to a pure dipole plus constant. 

In intermediate cases fits with both models (two 
Gaussians-|-quadrupole, Gaussian-|-dipole-|-quadrupole) , 
including a sextupole in the Fourier decomposition, pro- 
vide optimal descriptions and estimates of systematic un- 
certainty. Fits should be performed over the unbroken 2tt 
interval to respect and enforce periodicity of peak arrays. 
Analysis of SS or AS subintervals independently is incon- 
sistent with discrete Fourier series and can be misleading. 

VIII. SUMMARY 

I have examined the validity of jet structure iso- 
lated from dihadron azimuth correlations by combina- 
toric background subtraction based on ZYAM (zero yield 
at minimum) and conventional V2{pt) measurements. 
Simulations of ideal correlation data compared to pub- 
lished RHIC data demonstrate that in some cases ZYAM 
plus V2 subtraction produces substantial distortions and 
underestimation of jet yields. 

Combinatoric correlation analysis is a convenient 
method to access jet correlations in A-A collisions, but 
isolation of jet structure from the non-jet background 
presents technical challenges. The conventional back- 
ground model includes a constant offset Bq and an "ellip- 
tic flow" contribution nominally of the form u|(pt 1,^42) • 
The ZYAM principle is invoked to determine Bq, but 
ZYAM has no basis in conventional distribution analy- 



sis. The ZYAM method is defended by arguing that the 
detailed shapes of jet correlations are not known a pri- 
ori. Possible radiation effects and/or medium modifica- 
tions to fragmentation could distort jet-related structure. 
ZYAM is then said to reduce systematic errors arising 
from uncertainties about jet peak shapes. 

In this study I use simulations with known peak struc- 
tures to demonstrate that the ZYAM procedure can pro- 
duce large systematic errors. Systematic distortions may 
be much larger than any effects due to peak-shape un- 
certainties. Especially for more-central A-A collisions the 
ZYAM assumption leads to substantial underestimation 
of jet correlation amplitudes, interpreted in turn to imply 
large parton energy loss in a thermalized dense medium. 

I demonstrate that for typical correlations from RHIC 
collisions simple peak models adequately describe jet cor- 
relation data. Intrajet correlations (same-side peak) are 
well-described by a Gaussian. In special cases a modified 
Gaussian could accommodate width fluctuations due to 
radiation effects. Interjet correlations (awayside peak, 
back-to-back jets) are well-described by a Gaussian or, 
if the AS width is large (typical), by an azimuth dipole 
component. That model permits accurate determination 
of the background offset value, in contrast to ZYAM. 

Conventional measurements of V2 (pt) are susceptible to 
"nonffow" (jet correlations). In the conventional method 
of dihadron correlation analysis ^Kpt 1,^42) is overesti- 
mated, especially in more-central A-A collisions. As a 
result of f 2 oversubtraction the same-side jet peak is fur- 
ther reduced, and the away-side peak acquires a mini- 
mum at its center. Away-side distortions have been in- 
terpreted as manifestations of Mach shocks arising from 
sound propagation in a thermalized dense medium. 

As a result of ZYAM subtraction the impression may 
be formed that most partons are stopped in a thermal- 
ized dense medium ("opaque core"), that shock waves 
may be produced in the same medium as partons tra- 
verse the medium, giving rise to anomalous AS structure. 
In contrast untriggered angular correlations described by 
free model fits reveal that jet yields increase even faster 
than binary collisions in more-central Au-Au collisions. 
There is no isolated "away-side" jet. Away-side structure 
is undistorted and reflects jet-jet (parton-parton) corre- 
lations. Jets are modified in more-central collisions, but 
all partons survive as final-state correlation structure. 

This work was supported in part by the Office of Sci- 
ence of the U.S. DoE under grant DE-FG03-97ER41020. 
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